Sall genes regulate hindlimb initiation in mouse embryos.

Vertebrate limbs start to develop as paired protrusions from the lateral plate mesoderm at specific locations of the body with forelimb buds developing anteriorly and hindlimb buds posteriorly. During the initiation process, limb progenitor cells maintain active proliferation to form protrusions and start to express Fgf10, which triggers molecular processes for outgrowth and patterning. Although both processes occur in both types of limbs, forelimbs (Tbx5), and hindlimbs (Isl1) utilize distinct transcriptional systems to trigger their development. Here, we report that Sall1 and Sall4, zinc finger transcription factor genes, regulate hindlimb initiation in mouse embryos. Compared to the 100% frequency loss of hindlimb buds in TCre; Isl1 conditional knockouts, Hoxb6Cre; Isl1 conditional knockout causes a hypomorphic phenotype with only approximately 5% of mutants lacking the hindlimb. Our previous study of SALL4 ChIP-seq showed SALL4 enrichment in an Isl1 enhancer, suggesting that SALL4 acts upstream of Isl1. Removing 1 allele of Sall4 from the hypomorphic Hoxb6Cre; Isl1 mutant background caused loss of hindlimbs, but removing both alleles caused an even higher frequency of loss of hindlimbs, suggesting a genetic interaction between Sall4 and Isl1. Furthermore, TCre-mediated conditional double knockouts of Sall1 and Sall4 displayed a loss of expression of hindlimb progenitor markers (Isl1, Pitx1, Tbx4) and failed to develop hindlimbs, demonstrating functional redundancy between Sall1 and Sall4. Our data provides genetic evidence that Sall1 and Sall4 act as master regulators of hindlimb initiation.

Normal skeletal pattern formation in chick limb bud with a mesenchymal hole is mediated by adjustment of cellular properties along the anterior-posterior axis in the limb bud.

The chick limb bud has plasticity to reconstruct a normal skeletal pattern after a part of mesenchymal mass is excised to make a hole in its early stage of development. To understand the details of hole closure and re-establishment of normal limb axes to reconstruct a normal limb skeleton, we focused on cellular and molecular changes during hole repair and limb restoration. We excised a cube-shaped mass of mesenchymal cells from the medial region of chick hindlimb bud (stage 23) and observed the following morphogenesis. The hole had closed by 15 â€‹h after excision, followed by restoration of the limb bud morphology, and the cartilage pattern was largely restored by 48 â€‹h. Lineage analysis of the mesenchymal cells showed that cells at the anterior and posterior margins of the hole were adjoined at the hole closure site, whereas cells at the proximal and distal margins were not. To investigate cell polarity during hole repair, we analyzed intracellular positioning of the Golgi apparatus relative to the nuclei. We found that the Golgi apparatus tended to be directed toward the hole among cells at the anterior and posterior margins but not among cells at identical positions in normal limb buds or cells at the proximal and distal hole margins. In the manipulated limb buds, the frequency of cell proliferation was maintained compared with the control side. Tbx3 expression, which was usually restricted to anterior and posterior margins of the limb bud, was temporarily expanded medially and then reverted to a normal pattern as limb reconstruction proceeded, with Tbx3 negative cells reappearing in the medial regions of the limb buds. Thus, mesenchymal hole closure and limb reconstruction are mainly mediated by cells at the anterior and posterior hole margins. These results suggest that adjustment of cellular properties along the anteroposterior axis is crucial to restore limb damage and reconstruct normal skeletal patterns.

Cell-specific alterations in Pitx1 regulatory landscape activation caused by the loss of a single enhancer.

Developmental genes are frequently controlled by multiple enhancers sharing similar specificities. As a result, deletions of such regulatory elements have often failed to reveal their full function. Here, we use the Pitx1 testbed locus to characterize in detail the regulatory and cellular identity alterations following the deletion of one of its enhancers (Pen). By combining single cell transcriptomics and an in-embryo cell tracing approach, we observe an increased fraction of Pitx1 non/low-expressing cells and a decreased fraction of Pitx1 high-expressing cells. We find that the over-representation of Pitx1 non/low-expressing cells originates from a failure of the Pitx1 locus to coordinate enhancer activities and 3D chromatin changes. This locus mis-activation induces a localized heterochrony and a concurrent loss of irregular connective tissue, eventually leading to a clubfoot phenotype. This data suggests that, in some cases, redundant enhancers may be used to locally enforce a robust activation of their host regulatory landscapes.

SMAD4 target genes are part of a transcriptional network that integrates the response to BMP and SHH signaling during early limb bud patterning.

SMAD4 regulates gene expression in response to BMP and TGFß signal transduction, and is required for diverse morphogenetic processes, but its target genes have remained largely elusive. Here, we identify the SMAD4 target genes in mouse limb buds using an epitope-tagged Smad4 allele for ChIP-seq analysis in combination with transcription profiling. This analysis shows that SMAD4 predominantly mediates BMP signal transduction during early limb bud development. Unexpectedly, the expression of cholesterol biosynthesis enzymes is precociously downregulated and intracellular cholesterol levels are reduced in Smad4-deficient limb bud mesenchymal progenitors. Most importantly, our analysis reveals a predominant function of SMAD4 in upregulating target genes in the anterior limb bud mesenchyme. Analysis of differentially expressed genes shared between Smad4- and Shh-deficient limb buds corroborates this function of SMAD4 and also reveals the repressive effect of SMAD4 on posterior genes that are upregulated in response to SHH signaling. This analysis uncovers opposing trans-regulatory inputs from SHH- and SMAD4-mediated BMP signal transduction on anterior and posterior gene expression during the digit patterning and outgrowth in early limb buds.

Conserved and species-specific chromatin remodeling and regulatory dynamics during mouse and chicken limb bud development.

Chromatin remodeling and genomic alterations impact spatio-temporal regulation of gene expression, which is central to embryonic development. The analysis of mouse and chicken limb development provides important insights into the morphoregulatory mechanisms, however little is known about the regulatory differences underlying their morphological divergence. Here, we identify the underlying shared and species-specific epigenomic and genomic variations. In mouse forelimb buds, we observe striking synchrony between the temporal dynamics of chromatin accessibility and gene expression, while their divergence in chicken wing buds uncovers species-specific regulatory heterochrony. In silico mapping of transcription factor binding sites and computational footprinting establishes the developmental time-restricted transcription factor-DNA interactions. Finally, the construction of target gene networks for HAND2 and GLI3 transcriptional regulators reveals both conserved and species-specific interactions. Our analysis reveals the impact of genome evolution on the regulatory interactions orchestrating vertebrate limb bud morphogenesis and provides a molecular framework for comparative Evo-Devo studies.

Control of mouse limb initiation and antero-posterior patterning by Meis transcription factors.

Meis1 and Meis2 are homeodomain transcription factors that regulate organogenesis through cooperation with Hox proteins. Elimination of Meis genes after limb induction has shown their role in limb proximo-distal patterning; however, limb development in the complete absence of Meis function has not been studied. Here, we report that Meis1/2 inactivation in the lateral plate mesoderm of mouse embryos leads to limb agenesis. Meis and Tbx factors converge in this function, extensively co-binding with Tbx to genomic sites and co-regulating enhancers of Fgf10, a critical factor in limb initiation. Limbs with three deleted Meis alleles show proximal-specific skeletal hypoplasia and agenesis of posterior skeletal elements. This failure in posterior specification results from an early role of Meis factors in establishing the limb antero-posterior prepattern required for Shh activation. Our results demonstrate roles for Meis transcription factors in early limb development and identify their involvement in previously undescribed interaction networks that regulate organogenesis.

Generation of a new mouse line with conditionally activated signaling through the BMP receptor, ACVR1: A tool to characterize pleiotropic roles of BMP functions.

BMP signaling plays pleiotropic roles in various tissues during embryogenesis and after birth. We have previously generated a constitutively activated Acvr1(ca-Acvr1) transgenic mouse line (line L35) through pronuclei injection to investigate impacts of enhanced BMP signaling in a tissue specific manner. However, line L35 shows a restricted expression pattern of the transgene. Here, we generated another ca-Acvr1 transgenic line, line A11, using embryonic stem (ES) transgenesis. The generated line A11 shows distinctive phenotypes from line L35, along with very limited expression levels of the transgene. When the transgene is activated in the neural crest cells in a Cre-dependent manner, line A11 exhibits cleft palate and shorter jaws, while line L35 develops ectopic cartilages and highly hypomorphic facial structures. When activated in limb buds, line A11 develops organized but smaller limb skeletal structures, while line L35 forms disorganized limbs with little mineralization. Additionally, no heterotopic ossification (HO) is identified in line A11 when bred with NFATc1-Cre mice even after induction of tissue injury, which is an established protocol for HO for line L35. Therefore, the newly generated conditional ca-Acvr1 mouse line A11 provides an additional resource to dissect highly context dependent functions of BMP signaling in development and disease.

Developmental hourglass and heterochronic shifts in fin and limb development.

How genetic changes are linked to morphological novelties and developmental constraints remains elusive. Here, we investigate genetic apparatuses that distinguish fish fins from tetrapod limbs by analyzing transcriptomes and open-chromatin regions (OCRs). Specifically, we compared mouse forelimb buds with the pectoral fin buds of an elasmobranch, the brown-banded bamboo shark (Chiloscyllium punctatum). A transcriptomic comparison with an accurate orthology map revealed both a mass heterochrony and hourglass-shaped conservation of gene expression between fins and limbs. Furthermore, open-chromatin analysis suggested that access to conserved regulatory sequences is transiently increased during mid-stage limb development. During this stage, stage-specific and tissue-specific OCRs were also enriched. Together, early and late stages of fin/limb development are more permissive to mutations than middle stages, which may have contributed to major morphological changes during the fin-to-limb evolution. We hypothesize that the middle stages are constrained by regulatory complexity that results from dynamic and tissue-specific transcriptional controls.

Animals come in all shapes and sizes. This diversity arose through genetic mutations during evolution, but it is unclear exactly how these variations led to the formation of new shapes. There is increasing evidence to suggest that not all shapes are possible and that variability between animals is limited by a phenomenon known as "developmental constraint". These limitations direct parts of the body towards a specific shape as they develop in the embryo. Therefore, understanding the mechanisms underlying these developmental constraints could help explain how different body shapes evolved. The limbs of humans and other mammals evolved from the fins of fish, and this transition is often used to study the role developmental constraints play in evolution. This is an ideal model as there is already a detailed fossil record mapping this evolutionary event, and data pinpointing some of the genes involved in the development of limbs and fins. But this data is incomplete, and a full comparison between the genes activated in the fin and the limb during embryonic development had not been achieved. This is because most fish used for research have undergone recent genetic changes, making it hard to spot which genetic differences are linked to the evolution of the limb. To overcome this barrier, Onimaru et al. compared genetic data from the developing limbs of mice to the developing fins of the brown-banded bamboo shark, which evolves much slower than other fish. This revealed that although many genes commonly played a role in the development of the fin and the limb in the embryo, the activity of these shared genes was not the same. For example, genes that switched on in the late stages of limb development, switched off in the late stages of fin development. But in the middle of development, those differences were relatively small and both species activated very similar sets of genes. Many of these genes were pleiotropic, which means they have important roles in other tissues and therefore mutate less often. This suggests that the mid-stage of limb development is under the strongest level of constraint. Darwin's theory of natural selection explains that mutations drive evolution. But the theory cannot predict what kinds of new body shapes new mutations will produce. Understanding how the activity levels of different genes affect development could help to fill this knowledge gap. This has potential medical applications, for example, understanding why some genetic changes cause more serious problems than others. This work suggests that mutations in genes that are active during the mid-stage of limb development may have the most serious impact.

Exposure to tert-Butylphenyl Diphenyl Phosphate, an Organophosphate Ester Flame Retardant and Plasticizer, Alters Hedgehog Signaling in Murine Limb Bud Cultures.

Organophosphate esters have become widely used as flame retardants since the phase out of polybrominated diphenyl ethers. Previously, we demonstrated that some organophosphate esters, such as tert-butylphenyl diphenyl phosphate (BPDP), were more detrimental to endochondral ossification in murine limb bud cultures than one of the major polybrominated diphenyl ethers that they replaced, 2,2',4,4'-tetrabromodiphenyl ether. Here, we used a transcriptomic approach to elucidate the mechanism of action of BPDP in the developing limb. Limb buds collected from gestation day 13 CD1 mouse embryos were cultured for 3 or 24 h in the presence of vehicle, 1 µM, or 10 µM BPDP. RNA sequencing analyses revealed that exposure to 1 µM BPDP for 24 h increased the expression of 5 transcripts, including Ihh, and decreased 14 others, including Gli1, Ptch1, Ptch2, and other targets of Hedgehog (Hh) signaling. Pathway analysis predicted the inhibition of Hh signaling. Attenuation of Hh signaling activity began earlier and reached a greater magnitude after exposure to 10 µM BPDP. Because this pathway is part of the regulatory network governing endochondral ossification, we used a known Hh agonist, purmorphamine, to determine the contribution of Hh signaling inhibition to the negative impact of BPDP on endochondral ossification. Cotreatment of limbs with purmorphamine rescued the detrimental morphological changes in the cartilage template induced by BPDP exposure though it did not restore the expression of key transcription factors, Runx2 and Sp7, to control levels. These data highlight Hh signaling as a developmentally important pathway vulnerable to environmental chemical exposures.

IRX3/5 regulate mitotic chromatid segregation and limb bud shape.

Pattern formation is influenced by transcriptional regulation as well as by morphogenetic mechanisms that shape organ primordia, although factors that link these processes remain under-appreciated. Here we show that, apart from their established transcriptional roles in pattern formation, IRX3/5 help to shape the limb bud primordium by promoting the separation and intercalation of dividing mesodermal cells. Surprisingly, IRX3/5 are required for appropriate cell cycle progression and chromatid segregation during mitosis, possibly in a nontranscriptional manner. IRX3/5 associate with, promote the abundance of, and share overlapping functions with co-regulators of cell division such as the cohesin subunits SMC1, SMC3, NIPBL and CUX1. The findings imply that IRX3/5 coordinate early limb bud morphogenesis with skeletal pattern formation.

Transcriptomic similarities and differences between the limb bud AER and unique carapacial ridge of turtle embryos.

Evolutionary innovation may arise via major departures from an ancestral condition. Turtle shell morphogenesis depends on a unique structure known as the carapacial ridge (CR). This lateral tissue protrusion in turtle embryos exhibits similar properties as the apical ectodermal ridge (AER)-a well-known molecular signaling center involved in limb development. Still, how the CR influences shell morphogenesis is not entirely clear. The present study aimed to describe the CR transcriptome shortly before ribs were halted within its mesenchyme, as required for shell development. Analyses exposed that the mesenchymal marker VIM was one of the most highly co-expressed genes and numerous appendage formation genes were situated within the core of CR and AER co-expression networks. However, there were tissue-specific differences in the activity of these genes. For instance, WNT5A was most frequently assigned to appendage-related annotations of the CR network core, but not in the AER. Several homeobox transcription factors known to regulate limb bud patterning exhibited their highest expression levels in the AER, but were underexpressed in the CR. The results of this study corroborate that novel body plans often originate via alterations of pre-existing genetic networks. Altogether, this exploratory study enhances the groundwork for future experiments on the molecular underpinnings of turtle shell development and evolution.

HOX13-dependent chromatin accessibility underlies the transition towards the digit development program.

Hox genes encode transcription factors (TFs) that establish morphological diversity in the developing embryo. The similar DNA-binding motifs of the various HOX TFs contrast with the wide-range of HOX-dependent genetic programs. The influence of the chromatin context on HOX binding specificity remains elusive. Here, we used the developing limb as a model system to compare the binding specificity of HOXA13 and HOXD13 (HOX13 hereafter), which are required for digit formation, and HOXA11, involved in forearm/leg development. We find that upon ectopic expression in distal limb buds, HOXA11 binds sites normally HOX13-specific. Importantly, these sites are loci whose chromatin accessibility relies on HOX13. Moreover, we show that chromatin accessibility specific to the distal limb requires HOX13 function. Based on these results, we propose that HOX13 TFs pioneer the distal limb-specific chromatin accessibility landscape for the proper implementation of the distal limb developmental program.

The tissues and regulatory pattern of limb chondrogenesis.

Initial limb chondrogenesis offers the first differentiated tissues that resemble the mature skeletal anatomy. It is a developmental progression of three tissues. The limb begins with undifferentiated mesenchyme-1, some of which differentiates into condensations-2, and this tissue then transforms into cartilage-3. Each tissue is identified by physical characteristics of cell density, shape, and extracellular matrix composition. Tissue specific regimes of gene regulation underlie the diagnostic physical and chemical properties of these three tissues. These three tissue based regimes co-exist amid a background of other gene regulatory regimes within the same tissues and time-frame of limb development. The bio-molecular indicators of gene regulation reveal six identifiable patterns. Three of these patterns describe the unique bio-molecular indicators of each of the three tissues. A fourth pattern shares bio-molecular indicators between condensation and cartilage. Finally, a fifth pattern is composed of bio-molecular indicators that are found in undifferentiated mesenchyme prior to any condensation differentiation, then these bio-molecular indicators are upregulated in condensations and downregulated in undifferentiated mesenchyme. The undifferentiated mesenchyme that remains in between the condensations and cartilage, the interdigit, contains a unique set of bio-molecular indicators that exhibit dynamic behaviour during chondrogenesis and therefore argue for its own inclusion as a tissue in its own right and for more study into this process of differentiation.

Dynamic and self-regulatory interactions among gene regulatory networks control vertebrate limb bud morphogenesis.

Vertebrate limb bud outgrowth and patterning is controlled by two instructive signaling centers, the apical ectodermal ridge (AER) and the polarizing region in the posterior limb bud mesenchyme. Molecular analysis of limb bud development has identified a self-regulatory signaling system that operates between the AER and mesenchyme and orchestrates the dynamic progression of limb bud outgrowth and patterning. The first focus of this review are the gene regulatory networks (GRNs) and interactions that control the positioning of the fore- and hindlimb fields along the primary body axis, establish the initial axis polarity and control the precise positioning of the signaling centers. These early processes are largely controlled by activating and inhibiting interactions among types of transcriptional regulators expressed in specific territories. The second focus deals with the dynamic interactions among the GRNs that control limb bud patterning and outgrowth by responding to inputs from the self-regulatory limb bud signaling system. The final part describes the GRN interactions regulating digit morphogenesis and the Turing-type system that controls the periodicity of the digit ray pattern. This review highlights the significant progress made toward an integrative analysis and understanding of the morpho-regulatory systems that orchestrate patterning and outgrowth of vertebrate limb buds in time and space.

Gene Regulatory and Expression Differences between Mouse and Pig Limb Buds Provide Insights into the Evolutionary Emergence of Artiodactyl Traits.

Digit loss/reductions are evolutionary adaptations in cursorial mammals such as pigs. To gain mechanistic insight into these processes, we performed a comparative molecular analysis of limb development in mouse and pig embryos, which revealed a loss of anterior-posterior polarity during distal progression of pig limb bud development. These alterations in pig limb buds are paralleled by changes in the mesenchymal response to Sonic hedgehog (SHH) signaling, which is altered upstream of the reduction and loss of Fgf8 expression in the ectoderm that overlaps the reduced and vestigial digit rudiments of the pig handplate, respectively. Furthermore, genome-wide open chromatin profiling using equivalent developmental stages of mouse and pig limb buds reveals the functional divergence of about one-third of the regulatory genome. This study uncovers widespread alterations in the regulatory landscapes of genes essential for limb development that likely contributed to the morphological diversion of artiodactyl limbs from the pentadactyl archetype of tetrapods.

The Role of Retinoic Acid in Establishing the Early Limb Bud.

Retinoic acid (RA) was one of the first molecules in the modern era of experimental embryology to be shown capable of generating profound effects on limb development. In this review, we focus on the earliest events of limb development and specifically on the role of RA in establishing the domain of cells that will go on to form the limb itself. Although there is some consensus on the role of RA during the earliest stages of limb formation, some controversy remains on the mechanism of RA action and the requirement for RA signaling in forming the hindlimb buds.

Attenuated Fgf Signaling Underlies the Forelimb Heterochrony in the Emu Dromaius novaehollandiae.

Powered flight was fundamental to the establishment and radiation of birds. However, flight has been lost multiple times throughout avian evolution. Convergent losses of flight within the ratites (flightless paleognaths, including the emu and ostrich) often coincide with reduced wings. Although there is a wealth of anatomical knowledge for several ratites, the genetic mechanisms causing these changes remain debated. Here, we use a multidisciplinary approach employing embryological, genetic, and genomic techniques to interrogate the mechanisms underlying forelimb heterochrony in emu embryos. We show that the initiation of limb formation, an epithelial to mesenchymal transition (EMT) in the lateral plate mesoderm (LPM) and myoblast migration into the LPM, occur at equivalent stages in the emu and chick. However, the emu forelimb fails to subsequently proliferate. The unique emu forelimb expression of Nkx2.5, previously associated with diminished wing development, initiates after this stage (concomitant with myoblast migration into the LPM) and is therefore unlikely to cause this developmental delay. In contrast, RNA sequencing of limb tissue reveals significantly lower Fgf10 expression in the emu forelimb. Artificially increasing Fgf10 expression in the emu LPM induces ectodermal Fgf8 expression and a limb bud. Analyzing open chromatin reveals differentially active regulatory elements near Fgf10 and Sall-1 in the emu wing, and the Sall-1 enhancer activity is dependent on a likely Fgf-mediated Ets transcription factor-binding site. Taken together, our results suggest that regulatory changes result in lower expression of Fgf10 and a concomitant failure to express genes required for limb proliferation in the early emu wing bud.

Electrical stimulation of developmental forces reveals the mechanism of limb formation in vertebrate embryos.

Current knowledge on limbs development lacks a physical description of the forces leading to formation of the limbs precursors or "buds". Earlier stages of development are driven by large scale morphogenetic movements, such as dipolar vortical flows and mechanical buckling, pulled by rings of cells. It is a natural hypothesis that similar phenomena occur during limb formation. However it is difficult to experiment on the developmental forces, in such a complex dynamic system. Here, we report a physical study of hindlimb bud formation in the chicken embryo. We use electrical stimulation to enhance the physical forces present in the tissue, prior to limb bud formation. By triggering the physical forces in a rapid and amplified pattern, we reveal the mechanism of formation of the hindlimbs: the early presumptive embryonic territory is composed of a set of rings encased like Russian dolls. Each ring constricts in an excitable pattern of force, and the limb buds are generated by folding at a pre-existing boundary between two rings, forming the dorsal and ventral ectoderms. The amniotic sac buckles at another boundary. Physiologically, the actuator of the excitable force is the tail bud pushing posteriorly along the median axis. The developmental dynamics suggests how animals may evolve by modification of the magnitude of these forces, within a common broken symmetry. On a practical level, localized electrical stimulation of morphogenetic forces opens the way to in vivo electrical engineering of tissues.

Self-organized formation of developing appendages from murine pluripotent stem cells.

Limb development starts with the formation of limb buds (LBs), which consist of tissues from two different germ layers; the lateral plate mesoderm-derived mesenchyme and ectoderm-derived surface epithelium. Here, we report means for induction of an LB-like mesenchymal/epithelial complex tissues from murine pluripotent stem cells (PSCs) in vitro. The LB-like tissues selectively differentiate into forelimb- or hindlimb-type mesenchymes, depending on a concentration of retinoic acid. Comparative transcriptome analysis reveals that the LB-like tissues show similar gene expression pattern to that seen in LBs. We also show that manipulating BMP signaling enables us to induce a thickened epithelial structure similar to the apical ectodermal ridge. Finally, we demonstrate that the induced tissues can contribute to endogenous digit tissue after transplantation. This PSC technology offers a first step for creating an artificial limb bud in culture and might open the door to inducing other mesenchymal/epithelial complex tissues from PSCs.

Evolution of the avian digital pattern.

Variation in digit number has occurred multiple times in the history of archosaur evolution. The five digits of dinosaur limbs were reduced to three in bird forelimbs, and were further reduced in the vestigial forelimbs of the emu. Regulation of digit number has been investigated previously by examining genes involved in anterior-posterior patterning in forelimb buds among emu (Dromaius novaehollandiae), chicken (Gallus gallus) and zebra finch (Taeniopygia guttata). It was described that the expression of posterior genes are conserved among these three birds, whereas expression of anterior genes Gli3 and Alx4 varied significantly. Here we re-examined the expression pattern of Gli3 and Alx4 in the forelimb of emu, chicken and zebra finch. We found that Gli3 is expressed in the anterior region, although its range varied among species, and that the expression pattern of Alx4 in forelimb buds is broadly conserved in a stage-specific manner. We also found that the dynamic expression pattern of the BMP antagonist Gremlin1 (Grem1) in limb buds, which is critical for autopodial expansion, was consistent with the digital pattern of emu, chicken and zebra finch. Furthermore, in emu, variation among individuals was observed in the width of Grem1 expression in forelimb buds, as well as in the adult skeletal pattern. Our results support the view that the signalling system that regulates the dynamic expression of Grem1 in the limb bud contributes substantially to variations in avian digital patterns.